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Carbon monoxide is the third (after carbon dioxide and methane) most abun- 
dant carbon species in the atmosphere, but its current roles in the atmospheric and 
blogeochemlcal cycle of carbon are not quantitatively understood, Early work on car- 
bon monoxide included its discovery in the infrared solar absorption spectrum by 
Migeotte (1948) and an investigation its potential sources and sinks by Bates and With- 
erspoon (1652), The subject remained dormant for nearly two decades until the impor- 
tance of the hydroxyl radical in tropospheric chemistry was recognized (Levy, 1671; 
McConnell et al., 1971), The hydroxyl radical is derived from 

Og + hi/-* 0 e + 0( 1 D) (Rl) 


0( l D) + H b O -* ZV4 


(R7) 


the reaction 

CO + OH -♦ CO B + H (RIB) 

provides the only important sink for CO in the gas phase, and is responsible for deter- 
mining the chemical lifetime of atmospheric CO, In addition to destroying CO, the 
hydroxyl radical also destroys CH 4 

CH 4 + OH -» H b O + CH 3 , (R16) 

Reaction R16 initiates a sequence of reactions that ultimately lead to the production of 
CO, As first pointed out by McConnell et al, (1671), R16 represents a major source of 
CO, especially in the clean troposphere. It was argued (see Wofsy, 1676) that this 
source, together with that derived from incomplete combination of fossil fuels (Bates 
and Witherspoon, 1952), could account for the budget of atmospheric CO, However, 
Blermann et al. (1678) showed that reaction R18, which is primarily responsible for 
removing CO, is pressure-dependent in the presence of oxygen, and Bhould proceed 
nearly twice as fast as previously thought. A renewed search for sources of CO sug- 
gested that oxidation of non-methane hydrocarbons (NMHC), especially isoprenes and 
terpenes, emitted by trees (Zimmerman et al,, 1978) and biomass burning (Crutzen et 
al., 1679) could be important, However, it is difficult to quantify the magnitudes of 
these sources (see Logan et al„ 1681 for a detailed catalog and assessment of the 
uncertainties). 

The intense interest and effort on CO in the last decade results in the genera- 
tion of an extensive data base concerning its spatial distribution (see for example, 
Seiler, 1974) and seasonal variation (Dianov-Klokov et al., 1678). Previous studies of 
the CO cycle include photochemical models which neglected horizontal transport 
(Wofsy, 1976; Crutzen and Fishman, 1977; Logan et al., 1931), and Hameed and 
Stewart’s (1979) model, which considered north-south mixing by diffusion, However, 
the chemical lifetime of CO is believed to be at least several weeks, and hence both 
chemistry and transport should play a role in determining the CO concentration as a 
function of height, latitude and season. A study of the CO cycle with a three- 
dimensional global model, which explicitly accounts for atmospheric motions, should 
therefore be useful for helping to identify and quantify the sources and sinks of CO, 
The exercise is meaningful because an important part of our knowledge of CO is fairly 
well constrained; we can calculate the OH concentration in the troposphere to within a 
factor of 2 (OH determines the lifetime and methane oxidation source of CO), and we 
can estimate the anthropogenic (industrial) source with similar accuracy. The existing 
data sets offer roughly three independent pieces of information on CO; (a) latitudinal 
variation, (b) vertical variation, especially in the tropical upper troposphere and (c) 
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seasonal variation, at mid-latitude. In our model, theriJ is essentially only one free 
parameter, the magnitude of the plant (NMHC) source in the tropics, A large number 
of other parameters have also been tested, but they were shown to be unimportant, 
The concentrations of the hydroxyl radical are computed self-consistent ly. Since the 
lifetime of CO strongly affects its spatial and temporal variations, we consider our 
derived OH concentration as empirical values the* satisfy the CO observations, and 
their validity should transcend the assumptions of our simplified chemical model, 

In the following report, the essential results of our general circulation model 
study of CO are presented. 
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Abstract 


The carbon nonoxide cycle la studied by Incorporating the 
known CO sources and sinks In a tracer nodel which uses the winds 
generated by a general circulation nodel. The photochenlcal 
production and loss terns, which depend on OH radical 
concentrations, are calculated In an Interactive fashion. 
Comparison of the computed global distribution and seasonal 
variations of CO with observations is used to yield constraints 
on the distribution and nagnltude of the sources and s inks of CO, 
and the abundance of OH radicals In the troposphere. f Wc present' 
evidence that therfe"lB'"a substantial low latitude plant source of 1 
about 1.3 x 10*5 g yr“l, in addition to the anthropogenic and CH4 ) 
oxidation sourcesTj COur results suggest a globally averaged 
tropospheric OH concentration equal to 7 x 10 3 era -3 , and that 
this value might have been 20 percent higher In the 
pre-industrial atmosphere^ 


Introduction 


Carbon monoxide 1 s a frarr rnmpnnrnt ef the atmospher?; whose mixing 
ratio Is In the range 0.5 - 2.5 x 10~ 7 by volume (Migeotte, 1948; Seller, 
1974). The abundance and global distribution of the gas reflects Interaction^ 
between photochemistry, transport and surface sources and sinks. Quantltatlv^ 
knowledge of the CO cycle Is Important In Its own right, for the Information 
It Implies with regard to the global dispersal of atmospheric pollutants, and ( 
for an assessment of perturbations of the troposphere by man's activities. 

The identity and magnitudes of the major sources and sinks of CO are 
still uncertain. Suggestions for sources Include the Incomplete combustion of 
fossil fuels (Bates & Witherspoon, 1952), the oxidation of methane (McConnell 
et al., 1971), the oxidation of non-methane hydrocarbons (NMHC) emitted by 
forests (Zimmerman et al., 1978), marine biological activity (Seiler and 
Schmidt, 1975) and biomass burning (Crutzen et al., 1979). The primary loss 
mechanism for CO Involves the reaction (see Table la) 

CO + OH ♦ C0 2 + H (R18) 

where the hydroxyl radical f OH, Is derived from 
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• Oj ♦ hv ♦ 0 2 4 0<1 d> (II) 

o('W ♦ h k o -*» aoH 

(R7) 

(Uvy, 1971). Uptake by toll alcro-organlsas (Llabl and Sailer, 1976) may 
rapraaent an additional ainor alnk. 

Recent aatlaataa of the aean tropospheric OH abundance vary froa 1.5 * 

10® ea“ 3 (tfofsy, 1976) to 2 i 10® ca” 3 (Singh, 1977). Xaproved aatlaataa of 
tha OH concentration are needed to peralt aaeeasnent of the lapact of 
anthropogenic relaaaaa of halocarbona, auch aa CH*CC1« and CHjBr, on 
otratoapherlc 0*, aa veil aa for eat lasting the llfetlae of tropoepherlc 
pollutants, such aa SO* and nitrogen oxides (Yung et al., 1975; Sse, 1977, 
Haaeed at al., 1979). 

Previous studies of the CO budget and OH Include photocheaical aodela 
which neglected horlaontal transport (Vofsy, 1976; Crutzen and Flahnan, 1977; 
Logan at al., 1981), and Haaead and Stewart's (1979) aodel, which considered 
north-south nixing by diffusion. However, the chemical llfetlae of CO la 
believed to be at least several weeks, and hence both chemistry and transport 
should play a role In determining the CO concentration as a function of 
height, latitude and season. A study of the CO cycle with a three-dimensional 
global aodel, which explicitly accounts for ataospherlc notions, should 
therefore be valuable for helping to Identify sources and sinks of CO and the 
concentration of OH. 

Zn the following section we briefly describe the three-dimensional 
transport aodel used In our studies. In the subsequent section we describe 
aeveral experiments used to investigate the CO cycle. In the final two 
sections we discuss the results and Implications of the experiments. 

Description of the Model 

The numerical experiments of the carbon monoxide cycle were performed 
with a general circulation aodel (Hansen et al., 1981). The general 
circulation aodel is global In extent, with a horizontal resolution of 
8* x 10* (latitude x longitude). Seven layers, evenly spaced In 0 
co-ordinates, were employed between the surface and 10 ab, (31 km), with at 
Kjiast one layer entirely within the atratosphere at all locations. The aodel 
was run for five years, and the necessary dynamical variables for the fifth 
year stored off line for use In a "tracer aodel". Compared with 
climatological data, the three-dimensional wind fields, temperatures, and 
humidity, are all realistic. An Indication of the representativeness of the 
aodel 's wind fields Is presented In Pig. la for the surface winds In rebruary 
and In Pig. lb for the jet stream winds. The surface wind pattern shows the 
expected latitudinal variation, with easterlies at low latitudes, and cyclonic 
and antlcydonlc circulations at high latitudes. The southeast and northeast 
trade winds converge near the equator, the pattern being more apparent over 
tha ocean, where frictional affects are less than over land. Antlcydonlc 
patterns dominate over the oceans near 30* latitude, representative of the 
subtropic high pressure system. At about 50*N latitude the Aleutian low In 
the aid-Paclfic and the Icelandic low In the North Atlantic produce cyclonic 
wind patterns, while prevailing westerlies occur over the United States and 
Asia. In the aouthern hemisphere strong west winds are evident near 60*S. 

All these features predicted by our aodel are In accord with observations. At 
the jet stream level, strong winds are evident from 30* to 60*N, associated 
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with the subtropical and tha polar jata. Essential foaturaa include the 
•nhancad wind velocities near Japan and the east coaat of the United Staten, 
and a strong poleward coaponent over the far North Atlantic* At low latitudes 
ths jet stress is weak and predominantly easterly. Again, all such aodel 
predictions are consistent with observations. For purposes of this work, 
however, the aost relevant dynamical quantities are the meridional transport 
of water vapor, heat, aoaentua and angular aoaentua. A detailed coaparison of 
the values of the above quantities predicted by our GCM and the observations 
suggest that the agreeaent is better than 75Z. This provides strong evidence 
that the GCM can simulate the transport of trace species to an accuracy of 
75Z, adequate for the present investigation. At hourly time steps tha 
instantaneous aass distribution and the averaged horlsontal as as fluxes are 
fed into the tracer prograa. These aass fluxes are averaged over six hour 
intervals (Hahlaan and Moxla, 1978). The vertical aass fluxes are derived 
froa the horlsontal aass fluxes and the continuity equation. We also Included 
monthly-averaged vertical mixing by dry and aolst convection. Advectlon of a 
tracer is accomplished by a novel finite differencing technique. For each 
grid box wa predict four prognostic quantities: the aesn concentration of the , 

tracer n, and its spatial gradlant 9n. The echeme determines the four 
quantities for the next tine step using an upstream net hod. This algorithm is 
comparable in accuracy to a fourth order scheme (Mahlnan and Sinclair, 1977), 
but it has much less nuaerlcal noise. The nathenatical details are presented 
in Russell and Lerner (1981). 

The chemical sources and sinks of CO in the tracer continuity equation 
are provided by a simplified chemical model, which is, described in greater 
detail below. The concentration of the OH radical at each grid point is 
calculated using the current abundance of CO. We then use this value of OH to 
calculate the aethane oxidation source (R16) as well as the chemical sink 
(R18) for CO. The values of OH are updated monthly, which is sufficiently 
frequent, as indicated by comparison with test runs in which the updating is 
done every ten days. The non-photochemlcal sources and sinks are assumed to 
be time Independent, and modeled as constant release rates or destruction 
coefficients over the appropriate types of surfaces. 

Each experiment starts with a uniform distribution of CO, with a volume 
mixing ratio equal to 1 x 10”7. The tracer model is usually run for four 
years, repeatedly using one year's data generated by the general circulation 
model. That the tracer aodel has reched steady state in the fourth year is 
indicated by the lack of further change in the global CO dlstrubtlons between 
the third and fourth year. An a posteriori check on the mass balance in the 
tracer aodel reveals that it conserves mass to better than 1£ per year, which 
is adequate for the present investigations. 

Photochemical production and loss rates in the troposphere are 
constructed, using the set. of reactions shown in Table la. The tonally 
averaged CO distribution is used to generate photochemical equilibrium 
solutions for the free radicals OH, HO2, H2O2, 0( 1 D), NO and N02* 

Distributions of water vapor alxlng ratios and temperatures are taken from the 
general circulation aodel. Otone concent rat Iona and their seasonal variations 
are taken froa the observations of Herlng and Borden (1964). The 
concentrations of CH4 and H2 are set equal to 1.4 and 0*5 ppm, respectively, 
throughout the troposphere. The distribution of N0 X (NO + NO2) la taken froa 
Crutsen and Fishman (1977), with surface values of 0.1 ppb at 40*N and 
northward and 0.05 ppb southward, along with a acale height of 2 km. Because 
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of a lack of data, oo seasonal variations art prescribed for tha NO* 
coneant rations* Tha dlurnally-averaged phot odlsaoclat Ion rataa ara 
calculatad In tha usual way, with corractlons for Rayleigh tcattarlng (Yung at 
el*, I960)* Tha offset of eloud eovar la approximated, by dividing r.ll 
dissociation rataa in tha lower troposphere below the cloud heights (5 ka) by 
a factor of 2, alnce the average cloud cover la about 50Z. 

In arriving at the simplified chemistry shown In Table la we used several 
approximations regarding the Interactions between the hydroxyl radicals, 
nitrogen oxides and nethane-darlved radicala. Since we fix total NO, 

(■ NO + NOj), the only HOg - N0 X couplings are via the ratio, 

NO J 2 

wo 2 " 

and the formation of nitric add, 

N0 2 + OH 4 M ♦ RN0 3 + N . 

Since the primary fate of HN0 3 formed In (R12) Is loss by rainout, It Is a 
good approximation to assume that (R12) la a sink for HO^. The coupling 
between HO. and radicals derived from the oxidation of methane Is more 
complicated* We note that to an sxeellent approximation the destruction of 
each CH4 molecule leads to the production of one H 2 CQ molecule. Hence, the 
net result of methane chemistry on H0 X can be parameterized by the following 
expressions, 

P (H) - J4 IH 2 C0] - t 2 I OH] (2) 

P <H0 2 ) - k 20 [H 2 CO] {OH] + J« |H 2 CO] - Yl I OH] ^ + y 2 |OH] (3) 

L (OH) - k 20 IH 2 C0] |0H] - Yl l OH] ^ (4) 

whera 

Yl - *20 *16 t CH «l / + J 5 + *20 l° H l > 

Y2 " J 4 Yl / *20 

and P(x) and L(x) denote the production and loss rate (cm” 3 s _1 ) for species x 
respectively* 

Hlth these approximations we can calculate the photochemical equilibrium 
concentrations of the hydroxyl radicals (H, OH, HO* and H 2 0 2 ) by balancing the 
production and loss rates for these radicals, as given in Table la and 
expression (2)-(4>* This Involves solving four simultaneous equations. 
However, noting that 


( 1 ) 

(kl2> 


*18 tOOJ + *19 1*2^ + Y2 

t"l Hi t'6 2 r H I®»1 


( 5 ) 


and 
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i* g 


k 


!*2°2l 




J 3 4 k 15 (OB) + J 22 


|OH 2 ) 2 


(«> 


further simplifies the mathematics to seeking the solution for only two 
aon-llneer slgebrelc squetlons In two unknowns (OH snd HOj). The solution Is 
thsn used to luprove the set lusts of N0/N0 2 by (1). This procedure converges 
In s few lterstlons to better then 10"*. Tsble lb gives s represent stive 
sxsuple of the concent ret Ions of laportsnt redicel species celculeted by this 
slapllfled cheulcsl uodel. The spproxluete schene (l)-(6) csleulstes HQ* 
concent rst lone thet sre within 10X of those computed, using the exect schene 
given In Tsble Is. 

The strstoepherlc photochenlcel contributions ore token fron two- 
dlnentlonel nodel celeuletlons using the reectlon echene of Yung et si. 

(1980) , In which the lstltudlnel distribution of the long-lived species (Cl.. 
N?x. CH^. etc.) are prescribed by projecting the one-dlnenelonel nodel results 
along preferred nixing eurfecee. (HcElroy et al.. 1976). The distribution of 
stratospheric water vapor Is taken fron Wofsy (1976). 


Experlnents 

To investigate the CO sources and sinks, as well as the value of the neon 
tropospheric OH concent ret Ion. we performed eight expfcriments incorporating a 
wide range of production snd lose mechanisms. These experlnents are listed In 
Table 2. snd sre discussed below. 

In experlnent T1 we Included an anthropogenic source of 6.3 x 10** g/yr, 
(Seller. 1974), distributed over land areas in proportion to energy 
consumption. (Darastadter et al.. 1971). No seasonal variation In the 
anthropogenic source was Included as this effect Is thought to be small (Hall 
at al.. 1975). He also Included a source from nethane oxidation and a sink 
due to OH oxidation of CO. Both processes were distributed throughout the 
troposphere and stratosphere. 

Tracer experlnent T2 was the same as T1 except that the source fron 
nethane oxidation and loss by reaction with OH radicals were both arbitrarily 
reduced after each calculation to ensure a mean global abundance close to 
observed. This was done to determine the sensitivity of the nodel to the OH 
value. 

Experlnent T3 Included. In addition to the sources and sinks used in Tl. 
a plant source with a nagnltude of 4.3 x 10 1 * g/yr. equal to the nlnlnun 
astlnate of Zlnmernen et al. (1978). while their upper Unit of 1.3 x 10 1 ' 
g/yr, was used In T4. The latitudinal distribution of these sources was also 
taken fron Zlmmernan et al. for both experiments. Although the equatorial 
ground source has been paranetericed by using values estimated by Zlnmernen at 
al. (1978) for NMHC oxidation, there are other possibilities. For example, 
Crutxen at al. (1979) have suggested that the burning of vegetation either by 
forest fires or slash and burn agricultural practices, which are concentrated 
nalnly In the tropics, night be a large CO source. 
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Zn experiment T5, we used the anthropogenic source and included the 
■ethane oxidation eource and OH oxidation sink only in the stratosphere, to 
Investigate the Importance of transport to the stratosphere as a alnk, 

(Newell, 1977), and to laolate the affects of tropospheric photochealstry. 
Based on previous experience with this combination, we anticipated that an 
additional loss tarn would be necessary, thus, in experiment T6 we 
paraoeterlsed a aoll alnk by using a deposition velocity of 4 x 10*2 C n/sec, 
as estiaated by Liebl and Seller (1976). 

Zn experiaent T7, we Incorporated an ocean source, which is presumably 
associated with aarine biological activity. While estimates of a potential 
oceanic source vary widely, we have used the (conservative) value of 
0.4 x 10*4 g/yr, estimated by Seller and Schmidt (1975), which was distributed 
uniformly over the oceans. 

Finally, in experiment T8 we Included the same sources as in T4, except 
thst the anthropogenic source was excluded in order to Investigate the Impact 
of anthropogenic releases on the global distribution. 

The classification of sources into three groups- -anthropogenic, CH4 
oxidation and plant— is somewhat arbitrary and is done solely for convenience. 
For Instance, there would be a large anthropogenic component in the plant 
source, If a aajor component of that was due to slash and burn agricultural 
practices. 


Results and Discussion 

Figs. 2a and b depict vertically and zonally averaged CO concentrations 
as a function of latitude for the eight experiments, along with the observed 
values, (after Sellar, 1974). The experiment which produced the best fit was 
T4, which employed a maximum value of the plant source in addition to the 
anthropogenic and CH4 oxidation sources. All the other experiments resulted 
in values lower than observed except for T5, which, without any tropospheric 
chemistry, produced values twice as large as observed. Unless the transport 
to the stratosphere is substantially underestimated, this result indicates 
thst the stratosphere is not an Important sink. When the soil sink was added, 
as in T6, It reduced the concentration to a value less than observed. 

Obviously it would be possible to "tune" the deposition velocity by varying it 
within the uncertainty range of Liebl and Seller (1975), and therefore 
definitive conclusions about the strength of a soil sink can not be drawn from 
this experiment v*$,one. However, as will be discussed later, ‘the isotopic 
concentration of CO can be used to set constraints on the magnitude of the 
soil sink (Stevens et al., 1972). One can see that the experiments which lack 
any tropospheric chemistry, such as T5 and T6, or have reduced chemistry, such 
as T2, produce NH/SH ratios which differ from the observed by 202 or more. 

The qualitatively similar latitudinal distributions for T4 and T5 have 
different causes. Zn T4, the CO in the southern hemisphere is governed mainly 
by local sources and sinks (cf. Fig. Sa). In T5, CO rich air from northern 
mid latitudes is entrained into the rising branch of the Hadley circulation 
and is transported into the stratosphere, where it is destroyed. 

Stratospheric air which is poor in CO may then enter the mid-latitude southern 
hemisphere. Zn both cases, CO rich elr is transported southward and upward by 
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the northern hemispheric branch of the Hadley clrculat on end CO poor elr la 
transported northward and upward by the southern heals herlc branch, resulting 
in a sharp lnterhealspherie gradient et low altitudes. T8, which does not 
Include an anthropogenic source, has a poor lnterhealspherie difference. The 
latitudinal asyaaetry In T8 reflects the greater portion of land areas in the 
northern hemispheric tropics relative to the southern healaphere. T8 aay also 
describe the global CO distribution in the preindustrial atmosphere. The 
laportant assuaptlons and results of exji'iriments T1-T8 are sunmarlsed In Table 


The coaputed seasonal variations of CO at 55*N Is shown in Pigs. 3a and 
3b along with observations at Zvenlgorod, USSR (Dianov-Klokov et al., 1978). 
Pigs. 3c and 3d show the nodel derived seasonal variations at 20*N. The data 
is based on the observations of Seller et al. (1976) at Mauna Loa. Both data 
seta were shown to be relatively free of local pollution Influences, and 
should be representative of conditions in the clean troposphere. In general, 
all the models which had reduced OH levels In the troposphere could not 
satisfactorily reproduce the observed seasonal variation at 50*N, both in 
anplltude and phase. 

The observations of CO at 55*N and 20*N show little or no phase 
difference. This behavior Is In striking contrast to the seasonal variations 
of C 02 » which ahow significant differences between latitude belts, (Junge and 
Czeplak, 1968). The seasonal variations of CO 2 are driven aainly at high 
latitudes; at high northern latitudes the alnlaua occurs in late summer, which 
Is six months different from the southern hemispheric summer minimum. At low 
northern latitudes, with the Intrusion of air from the southern hemisphere, 
the phase Is deleyed. However, the observed seasonal variations of CO appear 
to be driven uniformly with latitude, Indicating a shorter lifetime and a 
dominance of local sources and sinks. This is consistent with seasonal 
variations In the OH concentration as the primary mechanism responsible for 
the seasonal variations of CO. The concentration of OH depends strongly on 
the duration and intensity of UV sunlight. Model derived OH values showed 
little or no variation In the times when maximum or minimum concentrations 
occur at different latitudes, the seasonal variations of CO observed by 
Seller et al. (1976) and Dianov-Klokov et al. (1978) likewise showed no 
difference In phase, a difference which would have been expected If the 
seasonal forcing were to arise from some feature such as a seasonal variation 
In the anthropogenic Input. 

If the seasonal CO data are* accurate this suggests that local chemistry 
aay be more laportant and the atmospheric OH may be even higher than It Is In 
the model. He have examined the sensitivity of the results In a model In 
which the OH concentrations were Increased by one third over the values In T4. 
The results showed that both the amplitudes and phases of the seasonal 
variations were in much better agreement with observations. Alternatively, 
there aay be a seasonal variation In the etrength of the low latitude source. 
As there Is little seasonality In primary productivity In tropics, which Is 
dominated by rain forests, we would expect no seasonal variation In the 
strength of the tropical source. If the scaling of hydrocarbon emissions to 
productivity used by Zimmerman et al. (1978) Is adopted. There aay be a 
seasonal variation in the CO source derived from vegetation burning. In the 
northern hemispheric tropics the dry season Is In winter (Richards, 1979). 
Enhanced forest fires In this season would contribute a source that could 
Improve the comparison between the model and observations. 




T5 and T6, which do not Include tropospheric chemistry, heve little 
seasonal variation* 12, with reduced chcalstry, hea a seasonal variation 
aoaewhat smaller than observed . The coaparlaona In Figs* 2-3 laply that both 
tropospheric chealatry and the anthropogenic source are necessary coaponenta 
of the CO budget. 

Ve also Investigated the latitude-height distribution of CO. Pig. 4 
ahovs an estlaate of the observed distribution, (Seller, 1974). The low level 
aaxiaua at upper aid-latitudes presumably due to the anthropogenic source, la 
apparent, as Is a lower latitude upper tropospheric bulge In the CO lsopleths. 
At 10*5, values Increase by sore than 752 between the surface and 12 tea. All 
the aodel experlaents fall to produce this feature except for T3 and T4, the 
experiments which Include the plant CO source. Pig* 4b shows the result froa 
year 4 for T2. The oalsslon of a low latitude eource Is prlaarily responsible 
for the absence of the equatorial bulge In this aodel. Fig. 4c presents the 
results for T4. The upper tropospheric bulge Is apparent as well as the 
aaxiaua at high latitudes. The distribution above the tropopause is 
unrealistic due to the failure of a seven-layer aodel to adequately reeolve 
the tropopause, and thus to accurately aodel troposphere-stratosphere 
exchange. 

It la of Interest to detera&ftty how the upper troposphere bulge Is 
produced In T4. As indicated In the above discussion, T4 concentrations are 
largely produced locally and this la true for the vertical distribution as 
well. The biologically produced CQ at low latitudes Is transported 
vertically, aalnly by the model's snail scale convection, in agreeaent with 
the aechsnlBB proposed by Falconer and Pratt (1980). Indication that a local 
ground source Is responsible for the observed bulge cones from the high 
altitude neasurenents of Gauntner et al. (1979) of CO and Altken nuclei over 
the tropics. Their results show that high CO concentrations are associated 
with high Altken nuclei counts, Indicating a common source for both of these 
quantities. The upper troposphere data of Seller (1974) were taken over the 
Atlantic, and are slallar to those observed by Gauntner et al. (1979) over 
Africa. In addition, the data of Gauntner et al. show higher values at 
equatorial latitudes over Africa than over the Pacific, The results of T4 are 
consistent with these observations. The longer distance froa the aajor 
northern hemispheric source regions result In lower values over t.he equatorial 
Pacific. In T2, the longer llfetlae of C allows for Its more complete 
redistribution around the globe. The Gauntner et al. data were obtained on 
only one flight, so aore data are needed to check If these are persistent 
features of the CO distribution, and to examine the seasonal, variability of 
the lnterhemlspheric transport of CO. 

Based on the complete set of experiments, we conclude that T4 Is the most 
successful combination of sources and sinks for simulating the observed CO 
distribution. A aore detailed analysis of 74 i»nd Its comparison with 12 will 
now be given. The continuity equation aay be written as: 


dM 

dt 


H+V+P-L+S 


(7) 


where N la Cha mi a of tracar, H the horizontal convergence associated with 
transport by wddles and tbs Man circulation, P and L the chaolcal production 
and loaa tana, V the vertical transport Into the troposphere froa the „ 
stratosphere, and 8 tha ground source, The results are shown In Pig. 5a and 
5b for T4 and T2, to emphasize tha roles of tha various teras in different 
experlMnts. T4 and T2 have soaewhat slallar latitudinal distributions as 
shown In Pig. 2a, pat tha additional source In T4, and tha reduced chemistry 
in T2 result in tha distributions being produced by different Mchanises. In 
T4 ouch of the CO is governed by a local cycle of sources and sinks, with the 
concentration peaking where the sources are a aaxlaua and in the Intervening 
30* latitude belt, where the horizontal convergence froa the neighboring zones 
is significant. In T2, the distribution at lower latitudes is no re strongly 
Influenced by the horizontal convergence. In this case, the long llfetiae of 
CO allows for its long range advection, but in so doing it also reduces auch 
of its variability. 

The results of these experlaents imply that the globally averaged OH 
concentration is approximately 7 x 10 5 cm -3 . (Singh (1977) has used 
observations of the lnterhemispherlc difference of Mthyl chloroform (CH3CCI3) 
to deduce globally averaged OH concentrations of at aost 4 x 10* cm - *. 

However, two recent etudles, Jeong and Kaufman, 1979; Kurylo et al., 1979), 
using Independent techniques, have derived a revised rate coefficient for the 
reaction CH3CCI3 + OH. This value is about a factor of two lower than that 
used by Singh in his analysis. These new results imply that his derived OH 
concentrations should be Increased by a factor of two and sire consistent with 
values obtained in experiment T4. Isopleths of the nean concentrations of OH 
in units of 10* molecules cm"* for the month of January and July computed in 
experiment T4 are shown in Figures 6a, b. These results for OH are in good 
agreement with those derived recently by Volz et al. (1981) on the basis of ** 
CO studies. 

Experiment T8 is designed to reconstruct the state of the atmosphere of 
the pre-industrial era. The input parameters for T8 are the same as those for 
T4 except for the industrial source, which is set equal to zero. A comparison 
between T4 and T8 suggests that there My have been a 652 change in the 
northern hemispheric CO concentration due to industrial pollution and a auch 
smaller change in the souther hemisphere. Indeed, a comparison between Shaw's 
(1958) data taken in Ohio during 1952-1953 and Dlanov-Klokov's data in the 
1970's in the USSR indicates that there might have been a substantial Increase 
in the CO concentration during this period. The mean value of OH in T8 is 
about 20X higher than that in T4, which suggests that there has been an 
Increase in the tropospheric lifetimes of species such as CH4, CH3CI , CH 3 Br 
and CHCI3. In our current investigation we have fixed the concentrations of 
CH4 and NO- at present atmospheric levels. If we had fixed the source 
strength or CH4, the shorter lifetime would have implied an Increase of CH4 
concentration in the atmosphere. Indeed, the recent observations by Rasmussen 
et al. (1981) provide strong evidence for tan Increase in atmospheric CH4. We 
have not conducted extensive studies on the sensitivity of our chemical model 
to Inputs such as N0 X and O3. The prlMry objective of this paper is to 
derive OH concentrations from the behaviour of CO in the model atmosphere. 

This study is, therefore, more like a diagnostic study of CO and OH than a 
prognostic one. Ve regard the observed distributions of O3 and N0 X as useful 
constraints in construction of chemical time constants for CO. It is 
precisely because of the uncertainties in the N0 n and O3 abundances in the 
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rat* coefficients for kay reaction, and In tha heterogeneous loss rataa that 
we hava adopt ad tha currant approach: a self-consistent atudy of the apatial 

dlatrlbutlon and taaporal variation* of CO and HO, to go beyond tha currant 
immediate obj active would be difficult to Justify. 

Global Maaureaenta of the '®0/ l8 0 ratio of CO by Stavana at al. (1972) 
typically show valuaa much lower than are derived from tha incomplete 
coabuetlon of foaall fuel*. Fractionation by destruction mechanisms cannot be 
invoked to explain tha data, Loat by reaction with OH radleala, R20, would 
favor tha lighter laotope (Wcinatock and Chang, 1974). Stavana at al. (1972) 
have measured tha fractionation in CO conauaption by aoll bacteria and found 
that it favora the lighter laotopaa. All of thnir northern healapherlc 
■aaauraaanta indicate a predominance of light leotoplc apeclea during the 
eumaar when aoll scavenging would be at ita maximum, which la contrary to what 
would be expected If the aoll were acting aa a major alnk. It would alao be 
difficult to explain the data by assuming equilibrium of CO with respect to 
uptake and release by the soil. By analogy with C0 2 , release froa the aoll 
and from plants may involve no fractionation, or again favor the lighter 
Isotope (Keeling, 1973). 

The models which do not Include vegetation sources have difficulty in 
explaining the laotople pattern found in the northern hemisphere, as they 
predict northern hemispheric CO to be mainly anthropogenic throughout the 
year. Exparlaaot T6, which assumes the only alnk to be aoll scavenging is 
also inconsistent with the Stevens et al. (1972) results aa the CO in the 
northern hemisphere would be even heavier than that released from combustion. 
As mentioned previously, the seasonal pattern of the > 8 0/ 16 0 ratio would be 
exactly opposite to what is observed. These results are consistent with the 
view that the soil plays a relatively minor role in the CO cycle. 

Since carbon-13 values from fossil fuel CO overlap those from natural 
sourcaa, these measurements may not be particularly useful as a diagnostic of 
source types. Measurements of carbon-14, on the other hand, might be more 
useful* The use of carbon-14 and oxygen-18 data would then minimise the 
uncertainties due to mixing of air masses. Figure 7 chows the annually 
averaged concentrations contributed by each of the sources we have included in 
experiment T4. They were derived by using each source separately with the OH 
field derived by T4 to determine the loses rate. The figure phows that 
virtually none of the anthropogenic CO emitted in the northern hemisphere 
reaches the southern hemisphere* The results also show that a alseable 
fraction of the CO at mid to high northern latitudes is not anthropogenic. 

Our results are in qualitative agreement with the major conclusions reached by 
Wcinatock and Niki (1972) and Stevens et al. (1972). More definitive 
statements cannot be made until more clean air measurements are available. 

The two sets of measurement* reported above were subject to local 
anthropogenic influences. 

Conclusions 

In this paper, we have examined an extensive set of experiments in the 
search for an understanding of four major features of atmospheric CO: (1) the 

meridional distribution, (2) the seasonal variation, (3) the equatorial bulge, 
and (4) the Isotopic composition. The results demonstrate the power and 
versatility of the general circulation model approach, as compared with more 
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traditional analyses (Wofsy. 1976; Crucian 6 Flshnan, 1977; Haaaad and 
Stewart. 1979)* The aost auceaaaful experiment. T4, which reproduces all the 
assantlal obsarvad features. requires. In addition to anthropogenic and CH 4 
oxidation aourcea, a large surface source of order 1.3 x I0 l * g/yr" 1 . This 
additional source oust span the low latitudes In the Hadley cell regime. In 
order to produce the observed equatorial bulge. 

Our aodel Is consistent with a globally averaged OH concent ret ion of 
about 7 x 10* cm"*. At this OH level, we are able to simulate the main 
features of the seasonal variation, observed by Dlanov-Klokov et al. (1978) In 
Zvenigorod (55*N). The seasonal data obtained by Seller at al. (1976) at 
Mauna Loa (20*N) laplles an even higher concentration of OH at lover latitudes 
than that In our T4 axperlaent or, alternatively, a possible seasonality in 
the plant aourca. The concentration of OH cannot, however, significantly 
exceed that calculated in T4» based on our current understanding of the 
CH 3 CCI 3 distribution (Singh. 1977; Chang and Fenner 1978; Jeong and Kaufman. 
1979; Kurylo at al., 1979). We are planning a general circulation aodel study 
of CH 3 CCI 3 with a 12-layer aodel in order to sake an Independent Investigation 
of tropospheric OH. 

Ho doubt a different and perhaps sore © iborate combination of sources 
and sinks and OH concentrations could also reproduce the observed features of 
CO. However, the aodel has been constructed with the fewest number of 
assumptions possible, and within the current small set of assumptions we are 
confident of the uniqueness of our results. We expect to laprove and refine 
our models as aore CO aeasurements become available.^ 

Measurements of the Isotopic content of the CO formed by various 
mechanisms such as vegetation burning, the oxidation of hydrocarbons such as 
Isoprenes and terpenes could be valuable In estimating the relative roles of 
tlvsve processes. Measurements of the seasonal variations of the total CO 
concentrations, and their Isotopic composition at remote locations could be 
used with data on the source characteristics to further quantify the 
magnitudes of the sources Involved, locations already In use for monitoring 
CO 2 , N 2 O and chlorofluorocarbons could be used for this purpose. Because the 
Intertroplcal convergence cone may be weaker ov.er the Pacific than the 
Atlantic (Atkinson and Sandler. 1970). It should present less of a barrier to 
lnterhemlspherlc exchange and local sources in determining seasonal cycles and 
interhemi spheric gradients of CO and other tracers in existing Pacific sites 
at the same latitudes. 
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Table la. Reactlone uaad for tha calculation of OH concentration in tha 

photehealcel aodal. Ttoo-and three -body rata coafflcianta 4. a given in 
unit* of co^a - * an d ca^e"*, raapactlvaly. Tha dlurnally-averaged aean 
dlasoclation ratal and ralnout ratal, (in unit! i” 1 ), ate given for 
aprlng equinox at 35*N, at tha ground. Tha value ahown for K12 la tha 
two-body high praaiura Halt. 



Reaction 

Rata Coefficient 

Reference 

R 1 

O 3 4 hv ♦ 0 2 4 0(*D) 

J1 - 1.92 x 10“ 6 , 

(a) 

12 

H0 2 4 hv ♦ NO 4 0 

J2 - 3.57 x 10 " 3 

(a) 

K3 

R 2 O 2 4 hv ♦ OH 4 OH 

J3 - 2.54 x 10 " 6 

<b) 

R4 

R 2 CO 4 hv ♦ HCO 4 H 

J4 - 5.04 x 10 “ 6 

(b) 

115 

h 2 co 4 hv ♦ h 2 4 CO 

J5 - 1.30 x 10-5 

(b) 

R 6 

0*0) 4 M ♦ 0( 3 P) 4 M 

k 6 - 4.0 x 10 " 11 

(b) 

R7 

0(*D) 4 R 2 0 ♦ OH 4 OH 

k7 - 2.3 x 10-*0 

(b> 

R8 

NO 4 O 3 ♦ N0 2 4 0 2 

k 8 - 2.1 x 10“22 ,-1450/T 

(a) 

R9 

OH 4 O 3 ♦ H0 2 4 0 2 

k9 - 1.5 x 10 T 12 .- 1000 /T 

(a) 

RIO 

H0 2 4 O 3 ♦ OH 4 20 2 

klO - 1.5 x ltf-l*a“ 600 ^ T 

(b) 

Rll 

H 4 0 2 + M ♦ H0 2 + M 

kll - 1.8 x 10"32 e 4340/T 

(c) 

R12 

OH 4 N0 2 ♦ HNO 3 

kl2 - 1.15 x 10 “ X1 

(a) 

R13 

OH 4 H0 2 ♦ H 2 0 4 0 2 

kl3 - 2.0 x 10 “ 10 

(d) 

R14 

H0 2 4 JP 2 ♦ H 2 O 2 4 O 2 

kl4 - 2.5 x 10 -12 

(b) 

R15 

H 2 0 2 4 OH ♦ H 2 0 4 H0 2 

kl5 - 1.0 x 10 -lla- 75 0 /T 

(b) 

R16 

CH 4 4 OH ♦ CH 3 4 H 2 0 

kl 6 - 2.3 x 10-12,-1710/T 

(b) 

R17 

NO 4 HO 2 ♦ OH 4 N0 2 

kl7 - 4.1 x 10 - 12 a 200 /T 

(b) 

R18 

C04 OH ♦ C0 2 4 R 

kl8 - 1.3 x 10-23(1 4 p ata ) 

( 0 ) 

R19 

H 2 4 OH ♦ H 2 0 4 H 

kl9 - 1.8 x 10"22 a -2400/T 

(b) 

R20 

R 2 CO 4 OH ♦ HCO 4 B 2 0 

k20 - 1.4 x 10 -11 

(c) 

R21 

HCO 4 O 2 ♦ CO 4 H0 2 

k21 - 5.7 x lO- 12 

(a) 

R22 

H 2 0 2 ♦ Ralnout 

J22 - 2.0 x 10*6 

(a) 


•NASA 1977. b NASA 1979. c Logan at al. (1978) 

d Hochanadal (1972) and DeMore (perional communication, 1978). 

•Wofiy, (1976). 




Table lb. A representative calculation of Lhe major special in the simplified 
chemical model for March at 35 # N at the ground. The kinetic rate 
coefficients are taken from Table 2a. All number densities are 
given In units of 


Input to Chemical Model 

Output from Chemical Model 

T - 281 *K 

[0 l D] - 1.66 x 10" 3 

M - 2.29 x 10 19 

(H] - 1.49 x 10” 1 

[H 2 0] - 2.38 x 10H 

[OH] - 5.93 x 10 5 

[ 0 3 ] - 8.4 x 10ll 

[H0 2 ] - 1.79 x 10® 

[NO x ] -1.55 x 10 9 

[H 2 0 2 ] - 3.49 x 10!0 

[CO] - 4.35 x 10 12 

[NO] - 3.61 x 10 s 
[N0 2 ] - 1.19 X 10 9 


Table 2. 

* 

Suaaary of sources and sinks used in the 
•odele. 

construction of the tracer 

Experiment 

t 

Sources 


Sinks 

T1 

Anthropogenic (a) 
CH^ Oxidation (b) 


CO 4- OH 

T2 

Anthropogenic (a) 

1/2 CH 4 Oxidation (b) 


1/2 (CO + OH) 

T3 

Anthropogenic (a) 

CH 4 Oxidation (b) 
Plant Source MIN (c) 


CO + OH 

T4 

Anthropogenic (a) 
CH 4 Oxidation (b) 
Plant Source Max (d) 


CO + OH 

T5 

Anthropogenic (a) 

CH 4 Oxidation Stratosphere 

(e) 

CO 4- OH Stratosphere (e) 

T6 

Anthropogenic (a) 

CH 4 Oxidation Stratosphere 

(e) 

CO 4- OH Stratosphere (e) 

T7 

Anthropogenic (a) 
CH 4 Oxidation (b) 
Ocean (g) 


CO 4- OH 

T8 

CH 4 Oxidation (b) 
Plant Source Max (d) 


CO 4- OH 

(a) 6.3 x 

10 1 * g/yr -1 (Seller. 1974) 




(b) based on OH calculated In the chealcal model (see text). 

(c) 4.3 x 10** g/yr -1 distributed according to Zlaoeraan et al. (1978). 
However we do not distinguish between NMHC and vegetation burning sources. 

(d) 1.3 x 10*3 g/yr -1 ,, otherwise sane as (c). 

(a) These reactlone take place in the stratosphere only. There Is no 
tropospheric chemistry. 

(f) deposition velocity over soil ■ 4 x 10 - * cas -1 . 

(g) 4 x 10 13 g/yr -1 (Seller 4 Schaldt. 1975). 
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Figure lb: Jet stream winds predicted by general 

circulation model for January. 
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Figures 3a- d: Calculated seasonal variations of CO for tracer experiments 

a) TI-T4 at 55 N, b) T5-T8 at 55 N observed values by, Dianov-Klokov 
et «1., 1978, c) T1-T4 at 20 N, d) T5-T8 at 20 N observed values by, 
Dianov-Klokov ct al.. 1978. 
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LATITUDE 

Figures 5a-b: Annual CO budget as a function of 

latitude, a) experiment T4, b) experiment T2. 

M - mass of tracer, H - horizontal flux covergence 
V - vertical transport lnot the toposphere from 
the stratosphere, P and L - photochemical production 
and loss terms and S - the oonchemlcal surface 
source terms. Units are lO 3 g/sec. 









experiment T4. a) Januarv, b) July. 
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